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We propose the sequential reaction process 15 0(p,7)(/3 + ) 16 as a new pathway to bypass of 
the ls O waiting point. This exotic reaction is found to have a surprisingly high cross section, 
approximately 10 10 times higher than the 15 0(p,/3 + ) 16 0. These cross sections were calculated after 
precise measurements of energies and widths of the proton-unbound 16 F low lying states, obtained 
using the H( 15 0,p) 15 reaction. The large (p, 7)(/3 + ) cross section can be understood to arise from 
the more efficient feeding of the low energy wing of the ground state resonance by the gamma decay. 
The implications of the new reaction in novae explosions and X-ray bursts are discussed. 
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Unbound nuclei play a major role in astrophysics. The 
proton-unbound 2 He and the alpha-unbound 8 Be nu- 
clei illustrate this fact. The former is involved in the 
p(p,P + )d reaction, first reaction of the pp chain of reac- 
tions governing the energy generation in the sun P, ■ 
The latter, whose lifetime is about 10 -16 seconds, is in- 
volved in the triple alpha reaction 0, which is at 
the origin of the formation of all the heavier elements. 
The proton-unbound nuclei 15 F and 16 F play an impor- 
tant role in X-ray bursts. These astronomical events arc 
known to happen in close binary systems, where accre- 
tion takes place from an extended companion star on the 
surface of a neutron star (type I X-ray burst). The ac- 
creted matter is compressed until it reaches sufficiently 
high pressure conditions to trigger a thermonuclear run- 
away, tn these explosive events, the carbon and nitrogen 
elements are mainly transformed into 1 O and 15 O by 
successive proton captures 0, || . Then, the pathway for 
new proton captures is hindered by the proton-unbound 
nuclei 15 F and 16 F. The reaction flux and the energy 
generation are then limited by the relatively slow (3 + - 
decay of 14 (t 1/2 =71 s) and 15 (t 1/2 =122 s), which 
create waiting points. The sudden and intense release of 
energy observed in X-ray bursts requires to circumvent 
the limited energy generation in breakout reactions. The 
15 0(cv,7) 19 Ne reaction is considered to be one of the key 
reactions in this context 0, • It makes the transition 
into the nucleosynthetic rp process (rapid proton cap- 
ture) which is responsible for an increased rate of energy 
generation and the synthesis of heavier elements In 
such explosive environments, 16 F is strongly populated in 
the ground state (g.s.) or in the first excited state, and 



leads to an equilibrium between formation and decay of 
this proton- unbound nucleus. From time to time before 
the proton is emitted, 16 F can capture anotherproton 
thus producing the 17 Ne particle stable isotope Q • This 
two-proton capture process was calculated to be signifi- 
cant for extreme densities (larger than 10 11 g/cm 3 ). In 
this letter, /3 + -decay of 16 F to 16 is proposed as an 
alternative channel. Two reactions 15 0(p,/? + ) 16 and 
15 0(p,7)(/3 + ) 16 are studied. Both reactions eventually 
proceed through the /3 + -decay of the intermediate un- 
bound 16 F g.s., which is fed directly by a proton cap- 
ture or indirectly through a proton capture to the first 
excited state followed by a 7-emission. The 7-decay oc- 
curs noticeably to the low energy wing of the g.s. reso- 
nance. Subsequent proton emission is dramatically hin- 
dered due to the fact that the low energy proton has to 
tunnel through the Coulomb potential of the 15 nucleus. 
These reaction channels have not been investigated so 
far and could speed-up the energy generation, competing 
with breakout reactions. The calculation of these reac- 
tion cross sections require the measurement of the ener- 
gies, widths, spins and parities of the low lying states of 
16 F. These were obtained from the measurement of the 
H( 15 0,p) 15 resonant elastic excitation function using 
low energy 15 O beam at the SPIRAL facility. 

The beam of radioactive 15 O nuclei was produced at 
the SPIRAL-GANIL facility through the projectile frag- 
mentation of a 95 A.MeV 16 primary beam on a thick 
carbon target. Mean intensities of 10 7 pps at an energy of 
1.2 A.MeV were obtained after post acceleration by the 
CIME cyclotron. A beam contamination of less than 1 % 
of 15 N was achieved using a vertical betatron oscillation 
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FIG. 1: Excitation function for the H( ls O,p) 15 reaction at 
180° in the center of mass frame. The line is a result of an 
R-matrix calculation using parameters from Table I. 



selection device ||| and a suitable degrader in the analy- 
sis line of LISE spectrometer [t| where the measurements 
were made. Two stable beams, 14 N and 15 N, were also 
used in similar experimental conditions for calibrations. 
The excitation function for the elastic scattering at these 
low energies can be described by the Rutherford scatter- 
ing, but shows "anomalies", i.e. various resonances that 
are related to individual states in the compound nucleus. 
The principle of the measurement is described in |lOl ITl| 
and references therein. A 31(1) /im thick polyethylene 
(CH 2 )„ target was used, thick enough to stop the beam 
inside. The scattered protons were detected by a silicon 
detector, placed at forward angles (180° in the center of 
mass frame) within an angular acceptance of 2°. Protons 
were identified using their energy and time-of- flight. The 
energy resolution was 3 keV in the center of mass frame. 
Fig. 1 shows the excitation function for the H( 15 0,p) 15 
reaction measured from 0.450 MeV to 1.1 MeV. The mea- 
sured cross section was reproduced by an R-matrix 
calculation using the code ANARKI which is seen 
to be in a good agreement with the data. A value of S p 
= -534 ± 5 keV was obtained for the proton separation 
energy in agreement with the recommended value [l4| . 
The R-matrix analysis was also used to extract the prop- 
erties of the first three states in 16 F, given in Table I. A 
significant difference was found between the present and 
the recommended value of the width for the first excited 
state [Tij . This width is an important parameter used in 
the calculations presented in the next section. 

The calculation of the 15 0(p,/3 + ) 16 cross section was 
made using the properties of 16 F g.s. resonance measured 
in the present work and the Breit-Wigner formula for a 
single- level resonance 0: 
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where A is the de Broglie wavelength, J are the spins, 
and En,TTot,^in,^out are the resonance energy, total 



TABLE I: Measured properties for the low-lying states in 16 F. 
E C m (keV) 6 E x (keV) n E x (keV) b JT T p (keV)" F p (keV) » 



534 ± 5 
732 ± 10 193 ± 6 
958 ± 2 424 ± 5 



0" 
198 ± 10 1" 
425 ±2 2" 



40 ± 20 25 ± 5 
< 40 70 ± 5 
40 ± 30 6 ± 3 



"Recommended values Il4l 
6 This work. 



width, and partial widths of the incoming and outgoing 
channels. In the (j?,/3 + ) case, Er — E g , s , the energy of 
the g.s. resonance, 
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the proton width, and 
T out corresponds to the /3 + -decay partial width. The 
energy dependance of the proton width T^ s - (E p ) for the 
incoming channel was taken into account by the using 
the relation: 



T^-(E p )=T^-(E g . s .) 



P(E P ) 
P(E g . s .) 



(2) 



where T^ s - (E g s ) is the proton width at the resonance 
energy and P(E p ) is the penetrability function under 
the Coulomb potential barrier. A partial lifetime for 
16 F((3 + ) of 1 second and a negligible branching ratio to 
the 15 0(p,/3 + ) 12 C+a final decay channel were assumed. 
This assumption is supported by the (3~ -decay properties 
measured in the mirror nucleus 16 N [3. The /3+-decay 
partial width was taken as a constant since the energy 
dependence of the Fermi function is small due to large 
Q (3 +=15417(8) keV [l4j. The calculated 15 0(p,/3+) 16 
cross section is shown in Fig. 2 as a function of the cen- 
ter of mass energy. The maximum of the cross section is 
observed at the energy of 534 keV corresponding to the 
16 F g.s. resonance. At this energy the (p, (3 + ) cross sec- 
tion is very small, about 10 -20 barns, since 16 F mainly 
decays by proton emission, which is ~ 10 20 times stronger 
than the /?+-dccay (since T^ s -{E g , s ) — 25 keV and r ou t 
= 0.66 10~ 18 keV). 




FIG. 2: 15 0(p,/3+) 16 and ls O(p, 7)(/3+) 16 cross sections 
are shown as a function of the center of mass energy. 

The calculation of the 15 0(p, 7)(/3 + ) 16 reaction was 
performed in a sequential manner, a schematic represen- 
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FIG. 3: Schematic representation of the ls O(p, 7)(/3 + ) 16 
reaction (see text). Two cases are represented. In red, 7- 
transitions populate the 16 F g.s. at the resonance energy. In 
that case 16 F mainly decays by proton emission. In blue, 
high energy 7-transitions populate the low energy wing of the 
g.s. resonance. In that case 16 F mainly decays by /3 + -ray 
emission. 



tation of this reaction is shown in Fig. 3. A proton cap- 
ture reaction to the first excited state of 16 F is considered, 
followed by a gamma decay to the g.s. resonance, from 
which a /3 + -decay branching ratio is taken into account. 
The cross section <j pl p(E p ) for the (p,j)((3 + ) reaction at 
the energy E p is an integration of the differential cross 
section over all possible energies of the gamma transition 
(since the g.s. has a large width): 



<J pl p(Ep) 



(Tpj(Ep,E 1 )P 1 (E 1 )Pp(Ep,E 1 )dE 1 (3) 



where a pl (E p , E 7 ) is the cross section to capture the pro- 
ton at the energy E p and to emit a 7-ray with an energy 
E 7 , P 1 {E 1 )dE 1 is the strength function 0,0], that is the 
probability for the 7-ray to have an energy between E 1 
and E~ ( +dE~ ( , and Pp(E p , E 7 ) is the branching ratio func- 
tion for the 16 F nucleus to decay by /3 + -ray emission. The 
first term a PJ (E p , Ej) is calculated using a Breit-Wigner 
formula with the following parameters Ei, Y\. AE p , E 7 ), 
T p (Ep), ri(J5 7 ) being the energy, total width, proton 
width and gamma width for the resonance correspond- 
ing to the first excited state of 16 F. The 7-ray is emitted 
from a 1~ state to the 0~ g.s., which corresponds to a 
Ml transition, whose energy dependence of the gamma 
width T\{E 7 ) is: 



T 1 JE 1 ) = T 1 JE 1 -E g . s .){ 
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(4) 



A gamma lifetime of 1 ps was obtained from the mir- 
ror nucleus 0], which corresponds to the partial width 
T^(Ei - E g . s .) = 0.6610-W. The strength function of 
the 16 F g.s. resonance was calculated assuming a Breit- 
Wigner parametrization: 



P 1 {E 1 )dE 1 



dE~, 



(5) 



and the normalization constant is: 

1 
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with AE = E v - E y 



E, 



g.s. 



and T Tot( E P 



T p (E p — E-y ) is the total width of the g.s. resonance. The 
j3 branching ratio is calculated using: 



P (E P1 E 7 ) 



Tp+r 9 p s -(E p -E y ) 



(7) 



Naively, one might have expected to obtain a small cross 
section for the (p, r y)(P + ) reaction, similar to the (p, (3 + ) 
one, since 7- and /3-widths are much smaller than proton- 
widths. Contrary to naive expectations, the (p, 7)(/3 + ) 
cross section is about 10 10 times larger than the (p, j3 + ) 
cross section, as shown in Fig. 2. The large ratio can be 
explained in the following way. As it has been shown pre- 
viously, there is only one (p,f3 + ) reaction for 10 20 (p,p) 
reactions. In the (p, "f){P + ) case, one 7-ray is emitted for 
10 8 incident protons (from the ratio of the widths) and 
about one 7-transition over 10 2 populates the low energy 
wing of the g.s. resonance (less than 50 keV above the 
proton emission threshold) where it is almost always fol- 
lowed by a /3 + -decay (P^ ~ 1). This implies that one 
incident proton over 10 10 induces a (p, j)(/3 + ) reaction, 
that is a factor 10 10 times larger than in the (p, (3 + ) re- 
action. 

In the following, uncertainties in the calculations and 
their evaluated effects on the results are discussed. The 
position and width of the low lying 16 F states were mea- 
sured with a high precision (see Table 1). The effect of 
the uncertainties in these measured parameters results 
in a change by less than a factor two in the calculated 
cross sections. The calculated (p, r y)(f3 + ) cross section is 
insensitive to the 16 F /3 + -decay lifetime, as a variation 
by a factor of 100 causes the cross section to change by 
only a factor of 2. The lifetime of the 7-transition is a 
sensitive parameter since the (p,7)(/3 + ) cross section is 
almost directly proportional to this parameter. A value 
measured in the mirror nucleus was used, but this as- 
sumption works only to within a factor of 10 [16]. The 
other excited states in 16 F were also studied and found to 
be negligible. The 15 0(p,7)(p,7) 17 Ne double indirect pro- 
ton capture reaction was not taken into account, neither 
the cross sections calculated, since it requires an appro- 
priate 3-body calculation. Moreover, cross sections may 
change by several effects, which remain to be evaluated 
as: non-resonant direct capture contributions, quantum 
interferences, continuum couplings [TtI ]. 

At a given temperature T of the gas inside the star, 
protons exhibit a Maxwellian distribution and the 
reaction rates < av > are calculated by integrating 
numerically the Maxwellian-averaged cross sections 
c(Ep) over all possible proton energies. The obtained 
reactions rates are shown in Fig. 4 (a) as a function 
of the temperature. The rate of the (p, (3 + ) reaction is 
negligible compared to that of the reaction (p,j){f3 + ) 
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FIG. 4: (a) 15 0(p,/3+) 16 and 15 0(p,7)(/?+) 16 reactions 
rates are shown as a function of the temperature. The 
15 0(a,7) 19 Ne reaction rate is also shown for compari- 
son, (b) Density versus temperature conditions where the 
15 0(p,7)(,3 + ) 16 reaction represents 10 to 50 % of the total 
reaction flux initiated by the 15 O nucleus. 



for all temperatures. To evaluate the impact of this 
latter reaction, it has to be compared with the com- 
peting /? + -decay of 15 and the 15 0(a,7) 19 Ne alpha 
capture reaction. Fig. 4 (b) shows the temperature 
and density conditions where the (p, 7)(/3 + ) reaction 
represents 10 to 50 % of the total reaction flux initiated 
by the 15 O nucleus. Boxes delimit conditions where 
novae and X-ray bursts can happen. For the lowest 
temperatures (< 10 8 K), the (p, 7)(/3 + ) reaction requires 
extreme densities (> 10 10 g cm -3 ) to compete with 
the 15 0(/3 + )-decay. For the highest temperatures (> 



1.1 10 9 K), the (a, 7) always dominates. In novae 
explosions, 15 nuclei mainly decay by /3 + -ray emission, 
the (p, 7)(/3 + ) reaction representing less than 1 % of the 
flux from 15 0. In X-ray bursts, the (p, 7)(/3 + ) reaction 
can represent up to 30 % of the total flux. Within the 
uncertainties of the calculations, the (p, j)(J3 + ) reaction 
could be faster than the (a, 7) reaction for temperatures 
up to 10 9 K. A more precise evaluation depends on the 
(a, 7) reaction rate (not well known) and on the relative 
abundances in hydrogen and helium, since one reaction 
consumes protons and the other alpha particles. In these 
extreme conditions, a new cycle of reactions is operating: 
15 0(p,7)(/3+) 16 0(p,7) 17 F(p, 7 ) 18 Ne(/3+) 18 F(p,a) 15 0. 
This new cycle could speed-up the CNO cycle and occur 
complementary to breakout reactions. The role of this 
new proposed cycle of reactions remains to be studied 
more carefully under various X-ray bursts conditions. 

In summary, it is shown that unbound nuclei can 
be involved in specific reactions that could play a role 
in astrophysics. Sequential (p, 7)(/3 + ) reaction, pro- 
ceeding trough an intermediate proton-unbound nu- 
cleus, was studied for the first time. The calculated 
15 0(p,7)(/3 + ) 16 cross section is found to be almost 10 10 
times larger than the direct 15 0(p,/3 + ) 16 reaction cross 
section. The large increase is mainly due to a strong 
feeding of the low energy wing of the 16 F g.s. reso- 
nance, where the subsequent /3 + -decay is favored. The 
{p-,l){P + ) could act in X-ray bursts, and would provide 
a steady burning scenario with a continuous depletion 
of 15 O. It is of great importance to study carefully the 
effects of this new reaction under various X-ray bursts 
conditions, and to demonstrate experimentally the ex- 
istence of (p,7)(/3 + ) reactions. The cross section of the 
15 0(p,7)(/3 + ) 16 reaction is calculated to be in the nano- 
barns range and can be measured using next generation 
intense RIB. More generally, several other unbound nu- 
clei as 19 Na or 15 F could also be involved in this type of 
reaction and remain to be studied. We thank the GANIL 
crew for delivering the 15 beam, M. Ploszajczak and A. 
Navin for stimulating discussions. This work has been 
supported by the IN2P3-IFIN-HH Program. 
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